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Abstract 

Deep-inelastic scattering events with a leading baryon have been detected by the HI ex- 
periment at HERA using a forward proton spectrometer and a forward neutron calorimeter. 
Semi-inclusive cross sections have been measured in the kinematic region 2 < Q 2 < 
50 GeV 2 , 6 • 1CT 5 < x < 6 • 1CT 3 and baryon p T < 200 MeV, for events with a final 
state proton with energy 580 < E' < 740 GeV, or a neutron with energy E' > 160 GeV. 
The measurements are used to test production models and factorization hypotheses. A 
Regge model of leading baryon production which consists of pion, pomeron and secondary 
reggeon exchanges gives an acceptable description of both semi-inclusive cross sections in 
the region 0.7 < E'/E p < 0.9, where E v is the proton beam energy. The leading neu- 
tron data are used to estimate for the first time the structure function of the pion at small 
Bjorken-x. 
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1 Introduction 



We report the measurement of the semi-inclusive cross sections for proton and neutron pro- 
duction in deep-inelastic scattering (DIS). The ep — > eNX data, where N represents either a 
final state proton or neutron, were obtained during 1995 and 1996 using the HERA accelerator 
at DESY where 27.5 GeV positrons collided with 820 GeV protons. Events containing lead- 
ing baryons were detected using the HI detector upgraded with a forward proton spectrometer 
(FPS) and a forward neutron calorimeter (FNC). 

The two multi-purpose experiments at HERA, HI and ZEUS, have observed a class of 
events which are characterized by the absence of final state particles in the region of phase 
space between the outgoing proton debris and the current jet [Jl],0]. These so called rapidity gap 
events contribute approximately 10% to the total DIS cross section and can be interpreted as 
being mainly due to interactions of a virtual photon with a pomeron In addition to virtual 

photon-pomeron interactions, one also expects meson exchanges to contribute to the total DIS 
cross section and to the production of leading protons and neutrons with small pr [f7|-[T0|]. Due 
to its small mass the most obvious candidate for such an exchange is the pion. If only pion 
exchange is responsible for leading baryon production and if isospin is conserved at the proton 
vertex, the ratio of neutron and proton production should be equal to two due to the difference 
in the Clebsch-Gordon coefficients for the n + n and ir p isospin | states. 

In this paper we present measurements of the semi-inclusive cross sections for proton and 
neutron production in the same kinematic region of x, Q 2 and pt, where x is the Bjorken scal- 
ing variable, Q 2 is the negative four-momentum squared of the virtual photon and p T is the 
transverse momentum of the final state baryon with respect to the beam direction. The mea- 
surements are compared to the results of a Regge model of leading baryon production and are 
used to test the pion exchange expectation for the ratio of neutron and proton production. We 
also compare our measurements to the predictions of the LEPTO 6.5 Monte Carlo program [ |Tl| ] 
which simulates baryon production without Regge dynamics by using an alternative formalism 
based upon soft colour interactions [0,0] and the string fragmentation model QTjn. Our cross 



section measurements are relevant for the determination of fracture functions which are a per- 
turbative QCD approach for describing the semi-inclusive production of hadrons in the proton 
fragmentation region [JT3|— |T7|] . 

We use the kinematic variables x, Q 2 and y to describe the inclusive DIS scattering process. 
They are defined as: 

~ q2 2 P-Q m 

x=- Q = -q y = r > (1) 

2p ■ q p ■ k 

where p, k and q are the four-momenta of the incident proton, the incident positron and the 
exchanged vector boson (7*) coupling to the positron. At ep centre-of-mass energy y/s they 
are related by Q 2 = sxy. 

The kinematic variables used to describe a final state baryon are: 

t = {p-p'f ~ - (1 -z){ ^ -nA z = \- q - {p - p>) ~ E'/E p , (2) 
z \ z p I q-p 
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where p' is the four-momentum of the final state baryon, m N is the mass of the final state baryon 
and m p is the proton mass. As defined, t corresponds to the squared four-momentum transferred 
between the incident proton and the final state baryon. It is different from the "rapidity gap" 
definition t — (p — py) 2 (used for example in reference [|]]), which defines the hadronic system 
Y at the proton vertex by the presence of a rapidity gap. For leading baryon production the 
definition of t in this paper is more appropriate even though it is only equal to the rapidity gap 
definition for events where the final state baryon is exclusively produced at the proton vertex. 

The four-fold differential cross section for baryon production can be parameterized by a 
semi-inclusive structure function, F^ 5 , defined by: 

dMep^eNX) _ W f, y> \ (3) 



dxdQ 2 dzdt xQ 4 \ tf 2[1 + R(x, Q 2 , z, t)] 

where a is the fine structure constant and R is the ratio between the absorption cross sections 
for longitudinally and transversely polarized virtual photons. In the kinematic range covered by 
this analysis, the structure function F^ B is rather insensitive to the value of R and we assume 
that R = 0. The variation from R = to R = oo leads to a 9% change at most in the resulting 
value of Fg ■* for the range of y covered in this measurement. 

The four-fold differential cross section integrated over < px < 200 MeV defines the 
semi-inclusive structure function F 2 v ' which we measure: 

d 3 a(ep — > eNX) f tmln Ana 2 ( y 



dx dQ 2 dz J t x Q 



(l-y + yW <4, (z,« 2 ,M)d« (4) 



where the integration limits are: 



J m% A (200 MeV) 2 
t min = -(l-z) M-mJ t = -- L +t mi n- (5) 

The structure function F^ 3 '' is denoted by F^f^ and F^ 3 - 1 for the semi-inclusive processes 
which have final state protons and neutrons respectively. We present measurements of F 2 in 
the range 0.73 < z < 0.88 and of F^^ for 0.3 < z < 0.9. In this paper we do not discuss the 
very high z region (z <: 0.95) which is most relevant for pomeron exchange and which has been 
used to measure the diffractive structure function F 2 ^ 3 ' ) [[3f-|(J]. 



2 HI Apparatus 



The central HI detector is described in detail in references Ql 8|— |2QQ . Here only the forward 
proton spectrometer, the forward neutron calorimeter and the crucial parts of the central HI 
detector used in this analysis are described. The coordinate system convention used by the HI 
experiment defines the positive 2-axis, or "forward" direction, as being that of the outgoing 
proton beam. The positive x-axis points towards the centre of the HERA ring. 
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The components of the central HI detector which are essential for this analysis are the 
backward electromagnetic calorimeter (SPACAL), the liquid argon calorimeter (LAr) and the 
central and forward tracking systems. The SPACAL calorimeter is used to determine the energy 
of the scattered positron whereas the hadronic final state is measured by the LAr calorimeter 
and the tracking systems. The SPACAL calorimeter has an electromagnetic energy resolution 



equal to a{E)j E = 7.1%/ y/ E[GeV]®l% as measured in an electron beam [ |19| ] and an angular 
acceptance of 153° < 9 < 177.5°. The LAr calorimeter has a hadronic energy resolution equal 



to a(E)/E ~ 0.5/ yE {E in GeV) as measured in a pion beam [ |21| ] and covers an angular 
range between 4° and 154°. Charged track momenta are measured using the central jet chamber 
(CJC), which consists of two concentric drift chambers covering an angular range between 15° 
and 165°, and by using the forward tracking system, which covers the angular range between 
7° and 25°. A uniform field of 1.15 T is produced using a superconducting solenoid which 
surrounds both the LAr calorimeter and the central tracking system. Luminosity is measured by 
detecting photons, from the Bethe-Heitler process ep — > ejry, in a crystal calorimeter situated 
at z = —103 m. 



2.1 Detection of Leading Protons 

The final state proton data were collected during 1995 with the HI forward proton spectrome- 
ter []2~2"|, [2~31] . In order to measure momenta of protons with scattering angles 6*^0.5 mrad with 
respect to the beam, the HERA machine magnets adjacent to the interaction region are em- 
ployed as spectrometer elements. Protons which have lost more than 10% of their energy in ep 
interactions appear after about 80 m at a distance of several millimeters from the nominal orbit 
so that they can be registered in detectors close to the circulating beam. The detector elements 
are mounted inside plunger vessels, so called Roman Pots, which are retracted during injection 
and are brought close to the beam after stable luminosity conditions are reached. The particles 
detected by the FPS are assumed to be protons. This is in agreement with Monte Carlo calcu- 
lations which show that approximately 98% of the charged particles observed in the FPS are 
protons. 

During 1995 HI operated two FPS stations, located 81 and 90 m away from the interaction 
point, which approach the beam from above. Each station is equipped with four planes of 
hodoscopes constructed of 1 mm scintillating fibres which are oriented ±45° with respect to 
the y axis. The scintillating fibre hodoscopes are 6 cm wide and 2.5 cm in height and are 
sandwiched between planes of trigger scintillators as sketched in Figure []]. Each scintillating 
fibre hodoscope has 240 scintillating fibres arranged in five densely packed layers. Four of the 
scintillating fibres, belonging to a common layer but separated by 11 fibres from each other, are 
attached to one cell of a 64 pixel position sensitive photomultiplier (Hamamatsu H-4139-20). 
Ambiguities in the hit combinations are resolved by using the segmented scintillators in front 
of and behind the hodoscopes. 

The FPS detectors measure space points with a resolution of a x = a y = 100 //m. There 
is no magnetic field between the FPS detectors and the space points measured at 81 and 90 m 
are used to define a track at 85 m. In both stations, ten out of the twenty layers have to show 
signals, in coincidence with trigger tiles, in order to reconstruct a track. The efficiency for a 
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minimum ionizing particle to produce a hit in a layer is typically 60% and the track reconstruc- 
tion efficiency is approximately 50%. With the help of the transfer matrices of the magnetic 
channel, which are known with high accuracy between the interaction point and 85 m, trajecto- 
ries are traced back to the interaction point and momenta are determined. Due to the fact that 
the HERA magnets lead to independent dispersion in both the horizontal and vertical planes, 
momenta can be measured twice by making use of the horizontal and vertical deflections. The 
two measurements have to agree within errors for tracks starting at the interaction point and this 
provides an efficient method for rejecting background tracks. 



2.2 Calibration of the Forward Proton Spectrometer 

In order to reconstruct momenta, the coordinates and slopes of scattered proton trajectories are 
transformed into a reference system with the beam direction as the primary axis. As the actual 
beam position is not known with the required accuracy of 0.5 mm, the actual beam orbit is 
determined for each fill using the nominal orbit as the first approximation. First the offset and 
tilt of the actual beam orbit with respect to the nominal one is determined in the horizontal 
plane by a fit that makes use of the fact that certain combinations of impact points and slopes 
at 85 m are "forbidden" for particles coming from the nominal interaction vertex [|22|]. For 
the calibration of the momentum measurement in the vertical plane, the difference between 
horizontal and vertical momentum measurements is used as additional input. Only protons, 
for which the momentum error arising from the uncertainty of the calibration in the horizontal 
plane is small, are used for this purpose. This procedure, which has been verified by Monte 
Carlo simulations of the FPS, leads to a momentum resolution of typically cr(E') = 6 GeV 
at 700 GeV. The angular resolution at the interaction point is cr(9 x ) = 5 /xrad, while <j(6 y ) 
varies between 5 and 100 /irad depending upon energy and angle. The calibration of the FPS is 
described in more detail in reference [(23]] . 



The FPS calibration was checked using high Q 2 DIS events with a forward rapidity gap [f2~2|]. 
For events with these characteristics the hadronic final state is well contained within the cen- 
tral HI detectors and one can compare the observed missing longitudinal momentum with the 
one measured by the FPS assuming that no particles escape in the forward region. The mean 
difference between the proton energy measured in the FPS and the energy expected from the 
calorimetric measurement is ( — 1 ± 9) GeV. From this we conclude that the scale of the FPS 
energy measurement is correct to within 10 GeV. 



2.3 Detection of Leading Neutrons 

During 1996 the HI experiment used a forward neutron calorimeter constructed of lead and 



scintillating fibres. The calorimeter, which was originally used by the WA89 experiment [24- 



260 at CEPvN, weighs approximately 10 tons and is located 107 m away from the nominal 
HI interaction point. Final state neutrons with production angles #<0.5 mrad are within the 
acceptance of the FNC. A schematic diagram of the FNC is shown in Figure |[ 

The forward neutron calorimeter consists of interleaved layers of 2 m long lead strips and 
scintillating fibres. The lead to fibre volume ratio is 4:1 and the nuclear interaction-length 
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Xj is 21 cm. The calorimeter is laterally segmented into hexagonal modules each of which is 
defined by coupling 1141 scintillating fibres to a common photomultiplier located at the rear of 
the detector. The height of a hexagonal module is 8.6 cm. A gap between the top and bottom 
parts of the calorimeter is necessary in order to have space for the proton beam pipe which 
passes through the calorimeter. There are 67 modules in the bottom part of the calorimeter and 
8 modules in the top part. 

The scintillating fibres are 1 mm in diameter and are orientated approximately parallel to the 
direction of the incident neutron. The attenuation length of the scintillating fibres is (1.7±0.2) m 
which has been measured using muons from cosmic events. Detailed GEANT [ ^7] ] simulation 
studies have shown that this small attenuation length is responsible for the high-energy tail 
(see Figure |J) that we observe in the neutron energy spectra. The high-energy tail is due 
to fluctuations in the longitudinal shower profile which lead to energy depositions close to the 
back-end of the calorimeter. Because of the small distance to the photomultipliers, the produced 
scintillation light is attenuated less than normal leading to an over-estimation of the incident 
particle's energy. The energy resolution of the calorimeter is a(E)/E ~ 20% for energies 
between 300 and 800 GeV. 

Two segmented planes of hodoscopes situated in front of the FNC are used to veto charged 
particles. Each plane is constructed of 1 cm thick hexagonal scintillator tiles which have the 
same lateral size as the calorimeter modules. The neutron detection efficiency of the FNC is 
(93 ± 5)%, the losses being due to coincidences in the veto hodoscopes which mostly originate 
from the back-scattering of charged particles produced during the neutron's hadronic shower. 
This efficiency was determined by measuring the rate of signals in the hodoscopes as a function 
of the radial distance away from the neutron impact position reconstructed in the calorimeter. 
Extrapolating the rate of signals in the hodoscopes to the region close to the impact position, 
the probability due to back-scattering was estimated. 

In this analysis we assume that all neutral clusters are produced by neutrons. Using the 
LEPTO Monte Carlo program t JTT| ] an d a GEANT [[27]] simulation of the HI beam line, we 



estimate that the background contribution due to other neutral particles (mostly photons and 
K®) is 6% for events with z > 0.2. For z > 0.6, the background contribution is 2%. 

All detector components, including the calorimeter and the hodoscope planes, are covered 
with lead sheets in order to shield them from synchrotron radiation. 



2.4 Calibration of the Forward Neutron Calorimeter 

The 75 modules of the FNC were initially calibrated at CERN using a 10 GeV incident electron 
beam. The FNC was positioned on a movable platform which allowed the response of each 
module to be measured separately. Preliminary calibration constants for the entire FNC were 
then determined by a matrix inversion procedure [|2"8"|]. After this initial calibration, the FNC 
had an approximately uniform response independent of impact position. 

After the calorimeter was installed in the HI beam line, run-dependent calibration con- 
stants were determined every few weeks by comparing the high-energy spectrum of neutrons, 
observed in interactions between the proton beam and residual gas in the beam pipe, with the 
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results of a pp — > nX Monte Carlo simulation based upon pion exchange [E^]. Hadronic 



pp — > nX data in the high z and low p T range, obtained at the ISR and by other experiments 
at CERN [^U|,pT]], are well described by pion exchange and have been used to constrain the 



pion flux factor [§|,f32|]. Since the pion flux factor determines the high-energy spectrum of final 
state neutrons, by comparing with the Monte Carlo simulation we are effectively calibrating the 
FNC with respect to previous experimental results. In the pp — ► nX Monte Carlo program, the 
acceptance and the energy response of the FNC are simulated by tracking neutrons through the 
GEANT Q simulation of the HI beam line. 

Figure |]a shows the uncorrected neutron energy spectrum observed in proton beam-gas 
interactions compared to the results of the Monte Carlo simulation. The two distributions are 
normalized to the same number of entries above 500 GeV. The peak position and the high- 
energy tail observed in the data are in good agreement with the simulation. Since the rate of 
neutron production with z < 0.5 is known to be underestimated by pion exchange [TJ2p, we 



do not use proton beam-gas data in this energy range for calibration purposes. This compari- 
son, between proton beam-gas interactions and the pion exchange Monte Carlo simulation, is 
the method we use to calibrate the FNC. We estimate a 5% energy scale uncertainty for the 
FNC based upon our comparison between proton beam-gas interactions and the pion exchange 
Monte Carlo simulation. 

Figure |3|b shows the same beam-gas energy spectrum compared to the neutron energy spec- 
trum observed in DIS interactions. Above 300 GeV the two distributions again agree very well 
in shape. This agreement supports the hypothesis that the pion flux factor is a universal property 
of the proton which is the same in both DIS and hadronic interactions Below 300 GeV, 

the sharp rise in the proton beam-gas energy spectrum is due to the trigger threshold used to 
obtain the data. 

The short-term gain variation of the FNC photomultipliers is measured by a LED monitor- 
ing system [ESI]. The light from seven light-emitting-diodes is coupled by optical fibres to the 



entrance windows of all the photomultipliers. The average response of the FNC photomultipli- 
ers to the LED light is used on a run-by-run basis to correct for small changes in the gain of the 
FNC photomultipliers. When there are stable beam conditions, the gain variation of the FNC 
photomultipliers is typically less than 0.1% during 30 minutes. 

The spatial resolution of the FNC was determined using charged particles and three small 
scintillator counters situated in front of the calorimeter. The scintillators are 3 x 3 x 10 mm 3 
and they are used in coincidence with a hodoscope tile to define a trigger. The spatial resolution 
of the FNC was determined to be: 

, N /5.13±0.81 , A 
^y(E) = ( ^fp^j + (°- 22 ± °- 07 )J cm ' ( 6 > 

where the reconstructed impact position was determined using the centre-of-gravity of the 
hadronic shower and an empirical formula which corrects for the hexagonal shape of the calorime- 
ter modules [^8L[33|]. 
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3 Event Selection and Data Analysis 



The final state proton and neutron data used to measure the semi-inclusive structure functions 
were collected during different years. The proton data were obtained during 1995 using a trigger 
which required a charged track through both detector stations of the FPS and a localized cluster 
in the backward (SPACAL) electromagnetic calorimeter. For part of the data a track candidate in 
the central jet chamber was also required by the trigger. During 1996, a trigger which required 
an energy deposit in the SPACAL electromagnetic calorimeter and the absence of out-of-time 
background signals was used to obtain the DIS data containing a high-energy neutron. 

During the offline analysis selection, criteria were applied to the data in order to reduce 
beam related backgrounds, events due to photoproduction and events from reactions in which 
the incoming positron lost a significant amount of energy by radiation. The DIS selection 
criteria used in the analysis are: 



A positron with energy E' e > 12 GeV in the angular range 156° <6 e < 177° was required 
which ensures that the scattered positron is within the acceptance region of the SPACAL 
electromagnetic calorimeter. 

The DIS kinematic variables were required to be in the range 2 < Q 2 < 50 GeV 2 , 0.02 < 
y < 0.6 and 6 -10 -5 < x < 6-10 -3 . The kinematic variables are reconstructed using the £ 



method [ p4| ] which combines the scattered positron energy and angle measurements with 
the quantity S, which is the sum over all hadronic final state particles of the differences 
between energy and longitudinal momentum. The E method has good resolution and 
keeps radiative corrections small over the entire kinematic range considered here. 

The quantity ~ Pz,i), which is calculated using the energy Ei and the longitudinal 

momentum p z ^ of all final state particles including the scattered positron, is expected to 
be twice the electron beam energy. This quantity was required to be > 40 GeV for the 
neutron analysis and > 41.6 GeV for the proton analysis. These cuts suppress radiative 
events and photoproduction background. 

The reconstructed vertex position was required to be within ±30 cm of the nominal vertex 
position in z. 

The proton analysis required at least one central track in the CJC with p T > 450 MeV to 
originate from the interaction vertex. 



Cuts related to the final state baryons are: 



• For the proton data, one forward track with 580 < E' p < 740 GeV and p T < 200 MeV 
was required to be detected by the FPS. Fiducial cuts on 9 X and 9 y , which depended on 
the proton energy, were applied to ensure that the track was observed in a region of the 
phase space where the acceptance was well understood and stable over the run period. 

• For the neutron data, one neutral cluster with E' n > 100 GeV and pt < 200 MeV was 
required to be reconstructed in the FNC. 
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After these cuts the data samples were grouped into 12 (x, Q 2 ) bins in the range 6 • 10~ 5 < 
x < 6 • 10~ 3 and 2 < Q 2 < 50 GeV 2 . The proton data sample consists of 1661 events and 
the neutron data consists of 10366 events. The total luminosities of the proton and neutron data 
samples are (1.44 ± 0.03) ptr 1 and (3.38 ± 0.07) ptr 1 respectively. 

The acceptances of the FPS and the FNC were determined by Monte Carlo programs in 
which protons or neutrons from DIS reactions were tracked through a simulation of the HERA 
beam line. The finite aperture of the beam line magnets limits the acceptance of both the FPS 
and the FNC. 

The FPS acceptance as a function of z was calculated using the RAPGAP Monte Carlo 



program which simulates pion exchange [p5|]. This Monte Carlo simulation gives a good de- 
scription of the shape of the uncorrected data as shown in Figure 0a for the observed proton 
energy spectrum. The FPS acceptance is approximately 80% for protons with 0.76 < z < 0.90 
and p T < 200 MeV. The LEPTO and ARIADNE @ Monte Carlo programs were used to 
check that the FPS acceptance is independent of the assumed production model and to estimate 
the systematic uncertainties. 

The corrected neutron energy spectrum was determined separately in each (x, Q 2 ) bin by 
using an unfolding procedure [B7p. The procedure uses Monte Carlo events to simultaneously 



correct the observed FNC energy spectrum for acceptance and migration effects. The LEPTO 
and RAPGAP Monte Carlo models were used to demonstrate that the unfolded neutron energy 
spectrum does not depend upon the production model used to correct the data. For neutrons 
with z > 0.4 and p T < 200 MeV the FNC acceptance is > 30%. 

Figures 0b and 0c show the observed neutron z and px spectra compared to the reweighted 
Monte Carlo simulation which results from the unfolding procedure. The data and Monte Carlo 
z distributions, shown in Figure 0b, are in agreement by construction since the unfolding pro- 
cedure does a fit to the data by reweighting the Monte Carlo as a function of z. The p T distribu- 
tions, shown in Figure 0c, demonstrate that the reweighted Monte Carlo gives a good description 
of a variable not used in the fit. 

The semi-inclusive structure functions F^f^ and F^ 3 ^ have been corrected to the Born 



level. Radiative corrections were calculated using the program HERACLES [|38|]. In all (x, Q 2 ) 
bins they are smaller than 6%. We have included a 2% systematic error on our radiative cor- 
rections due to hadronic corrections and higher order processes which are not simulated by the 
HERACLES code. 

There are three types of systematic errors: normalization errors, errors which depend on the 
final state baryons and errors which are different for each (x, Q 2 ) bin: 



• Normalization Systematic Errors 

- For the proton analysis the normalization error is 5.6%. The main contribution to 
this error is the 5.0% uncertainty in the proton reconstruction efficiency. 

- The normalization error for the neutron analysis is 5.7%. The 5.4% uncertainty in 
the neutron detection efficiency is the largest component of this error. 
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The systematic errors on the total integrated luminosities, which are approximately 2%, 
are included in the normalization uncertainties. 

• Final State Baryon Systematic Errors 

- For the proton analysis, these errors are between 4.8% and 19%. Errors in the mi- 
gration corrections for the proton energy intervals, which depend on the accuracy of 
the calibration procedure, were evaluated from Monte Carlo studies and range be- 
tween 4.5% and 19% for the different bins. Additional errors due to the uncertainty 
in the acceptance of the the fiducial cuts were evaluated by comparing the results 
obtained using the RAPGAP, LEPTO and ARIADNE Monte Carlo models. 

- The FNC energy spectrum, after being corrected for acceptance and migration ef- 
fects by using the unfolding procedure, has large systematic errors. We have varied 
the FNC energy scale by ±5% and reweighted the Monte Carlo data as a function 
of pt to estimate these systematic errors. We have also used the shape of the im- 
pact point distribution and its maximum, which defines the zero degree direction, 
to determine the FNC acceptance. A systematic error is applied corresponding to 
the difference between this acceptance method and the acceptance determined using 
the unfolding procedure. These uncertainties lead to systematic errors which range 
between 16% and 58% for the corrected neutron energy spectrum. 

• Errors which depend on x and Q 2 

- The systematic errors which differ in each (x, Q 2 ) bin range between 7.7% and 13% 
for the proton analysis. The uncertainties in the trigger efficiency and in the correc- 
tions for migrations between different bins in x and Q 2 are the main contributions to 
these systematic errors. The uncertainty in the acceptance and in the migrations as 
a function of x and Q 2 was determined by varying the energy scale of the SPACAL 
calorimeter by ±1.5%, by varying the electron scattering angle by ±1 mrad and 
by simulating, according to our knowledge of the hadronic energy scale of the LAr 
calorimeter, a 4% uncertainty in the measurement of E. 

- For the neutron analysis, the systematic errors which depend upon x and Q 2 range 
between 3.1% and 7.2%, since the measurement of the scattered positron was further 
improved in 1996. The main source of these systematic errors is the uncertainty in 
the acceptance and in the migrations as a function of x and Q 2 which was determined 
by varying the energy scale of the SPACAL calorimeter by between ±1% and ±3% 
(depending upon the energy of the scattered positron), by varying the reconstructed 
angle of the scattered positron by ±0.5 mrad and by varying the energy scale of the 
LAr calorimeter by ±4%. 

These errors can be compared with the statistical ones. The statistical errors for the proton 
analysis lie between 9.6% and 30% for 90% of the data points. For the neutron analysis, the 
statistical errors range between 4.7% and 29% for 90% of the data points. 
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4 The Semi-Inclusive Structure Functions F^ 3 ^ and F^ 3 ^ 

Our measurement of the semi-inclusive structure functions F^ 3 ^ and F^ 3 '', for leading pro- 
tons and neutrons with pr < 200 MeV, are shown in Figures |5| and ^]. The inner error bars show 
the statistical errors and the full error bars show the statistical and systematic errors added in 
quadrature. Tables [1] and ^ list the values of the semi-inclusive structure functions shown in the 
figures. The data are compared to the predictions of the LEPTO and RAPGAP Monte Carlo 
models (TJ§. 

The LEPTO 6.5 Monte Carlo program simulates baryon production using soft colour in- 
teractions and the JETSET string fragmentation model jJTl|,[T4t]. Soft colour interaction mod- 
els have been proposed to explain large rapidity gap events and the production of final state 
baryons [JT^,[13j] . In these models, the colour structure of the partons interacting with the virtual 
photon is modified by non-perturbative soft gluon exchanges which can lead to the production 
of colour neutral partonic subsystems separated in rapidity. After the fragmentation process, a 
high-energy baryon may be produced separated by a large rapidity gap from the remainder of 
the hadronic final state. 

The LEPTO Monte Carlo model describes the general shape and magnitude of the neutron 
data over the entire z range. It fails however to describe the rate of leading proton production 
and the rise in F^ 3 '' as a function of Q 2 . The leading order parton distributions for the proton 



by Gliick, Reya and Vogt (GRV) [g^gg], and the default value of 0.5 as the probability for soft 
colour interactions, were used to calculate the LEPTO Monte Carlo predictions. 

The RAPGAP Monte Carlo program simulates leading baryon production using pion ex- 
change. In the Monte Carlo program, the cross sections for leading proton and neutron produc- 
tion are proportional to the product of the pion flux factor and the pion structure function. The 
pion flux factor determines the energy and p T spectra of the final state baryons and is identical 
for proton and neutron production except for a factor of two. For the RAPGAP Monte Carlo 
predictions shown in Figures |5| and || we have used the pion flux factor determined by Holtmann 
et ah [|8|]. The rate of leading baryon production depends also upon the values of the pion parton 
distributions and we have used the leading order parametrization by GRV [BPLHO]. 



The RAPGAP Monte Carlo gives a reasonable description of the high-energy neutron data 
with z > 0.7 but it fails to reproduce the absolute rate of proton production. In the low-energy 
region where the final state neutron has < 70% of the incident proton's energy, the RAPGAP 
Monte Carlo program is not valid since additional physical processes, not simulated by the 
program, are expected to contribute significantly to the production of neutrons. 

It is interesting to compare the magnitudes of F^ 3 '' and F^ 3 ^. For z > 0.7, the semi- 
inclusive cross section for proton production is larger than the cross section for neutron pro- 
duction in any specific (x, Q 2 ) bin. This result rules out pion exchange as the main production 
mechanism for leading protons since pion exchange models predict that the ratio of neutron and 
proton production should be equal to two. 
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5 Factorization and Scaling Violations of F^ 3 ^ and F^ 3 ^ 

Presuming that leading baryons emerge from reactions where the virtual photon is absorbed by 

LB (3) 

a colourless object inside the proton, the structure function F 2 v ; should factorize into a flux 
factor f(z) which is only a function of z, and a structure function F^ 2 ^ which depends upon 
[3 and Q 2 . The quantity (3 = x/(l — z) may be interpreted as the fraction of the exchanged 
object's momentum carried by the quark or gluon interacting with the virtual photon. 

Alternatively one may assume models [ |T2| ] which are not based on the exchange of colour- 
less objects so that (3 can no longer be interpreted as a momentum fraction. In such scenarios, 
one might expect factorization in the variables x, Q 2 and z, if the deep-inelastic scattering 
process off the proton is independent of the proton fragmentation. The "hypothesis of limiting 
fragmentation" which states that target fragmentation is independent of the incident pro- 
jectile's energy, also implies that final state baryons emerge from a process which is insensitive 
to x and Q 2 . 

To test both factorization hypotheses, fits were made to the proton and neutron data sepa- 
rately assuming the following general forms for F^ 3 ^ [|2"3|,|2"9"[]: 

F^ 3) (f3,Q 2 ,z) = f(z)-F^ 2 \f3,Q 2 ) (7) 

Ff (3) (x, Q 2 ,z) = f(z)-F^ m (x, Q 2 ), (8) 

where the discrete-function f(z) is expressed by three free-parameters and for F^ 6 ^ Q 2 ) in 
equation a functional form, based on the leading terms of a phenomenological parametrization 
of the proton structure function [03], was chosen. For F 2 B ^(:r, Q 2 ) in equation |S|, (3 was 
replaced by x. 

The data are consistent with both factorization hypotheses and the fit results yield simi- 
lar x 2 /ndf. A possible explanation for this result is that F^ 2 - 1 is proportional to the proton 
structure function which for x < 0.1 is of the form F 2 ~ x~ x ( Q ) p2[]. Since (3 and x are 
highly correlated and have similar magnitude due to the restricted range of z, this also implies 
that F^ 2 ^ ~ (3~ x( - q2 \ The data have therefore relatively limited sensitivity to a difference of 
factorization in these two variables. 

In order to quantify the scaling violations observed in the data, F^f^ and F^ 3 ^ have been 
fitted separately for each fixed value of (3 to the form: 

F^ 2 \(3,Q 2 ) = a((3) + b(P)-\ogQ 2 , (9) 

with Q 2 in GeV 2 . The values of b((3) /F 2 B(2) , which are a measure of the scaling violations, are 
plotted in Figure ^. Only the fit results which arise from the proton data with z = 0.732 and the 
neutron data with z = 0.7 are shown since the results from the other z values are similar. The 
measurements of F 2 P( - 3 ^ and F 2 N ^ in the lowest x bin have not been used since there is only a 
single Q 2 value. The data are compared to the scaling violations dF 2 /d(logQ 2 )/F 2 predicted 
and observed in the inclusive structure functions of the pion and proton respectively. The pion 
and proton structure functions have been calculated using the GRV leading order parametriza- 
tions P^,HUp. The scaling violations observed in the semi-inclusive structure functions F 2 



and F 2 are similar in size and shape and are close to those seen in the GRV parametrizations 
of the pion and proton structure functions. 



14 



6 Comparison to a Regge Model of Baryon Production 



Assuming a simple Regge expansion and the dominance of a single Regge exchange, the dif- 
ferential cross section for leading baryon production as a function of z at fixed t should be 
proportional to (1 — z)~ n . Here n = 2a(t) — 1, and a(t) specifies the Regge trajectory of the 
dominant exchange. For the leading neutron data with 0.7 < z < 0.9 shown in Figure ^ the 
falling z spectra suggest a value of n, averaged over the t-dependence of the baryon production 
cross section, which is approximately equal to —1. This implies that the average value of a(t) is 
consistent with zero which is naively the expectation of pion exchange. In contrast the leading 
proton data, shown in Figure |5], do not depend strongly on z so that the average value of a(t) is 
larger than the value suggested by the neutron data. This is consistent with the dominance of a 
trajectory with the intercept a(0) ~ 0.5 which was found to be the sub-leading contribution in 
the diffractive region at larger z [|5p. 

Figure |] shows a comparison between the leading baryon structure functions with 0.7 < 
z < 0.9 and a Regge model of baryon production. In the model, the contribution of a specific 
exchange i is determined by the product of its particle flux fi/ p (z, t) and its structure function 
F£ evaluated at (j3, Q 2 ). For leading baryon production with p T < 200 MeV we therefore have: 

F™ {3 \p,Q 2 ,z) =J2 (j t tmm fi/ P (z,t)dt\ -F 2 (/3,Q 2 ), (10) 

where i denotes the pion, the pomeron and secondary reggeons (for example p, u>, a 2 and f 2 ). 
The integration limits t and t min are given by equation |5J. 

In the Regge model, we assume that the neutral pion, the pomeron and the f 2 all contribute 
to leading proton production. We neglect the contributions due to the other secondary reggeons 
because there is no sensitivity to them in the data, and because they have been estimated to be 



much smaller than the contribution due to f 2 exchange p3| , |44| ]. A comparison of total hadronic 
cross section measurements has resulted in the estimate that the flux of reggeons which have 
isospin equal to one (p and a 2 ) is only « 3% of the flux of reggeons with isospin equal to zero (u 
and f 2 ) [ff3]]. Regge phenomenology also predicts f 2 dominance, among isoscalar trajectories 
in the present case, in contrast to exchange degeneracy for elastic scattering processes [fflp. 

For leading neutrons, we assume that they are produced by charged pion exchange only. 
In the limited pt range of the data, leading neutron production due to p and a 2 exchanges 
has been estimated to be more than an order of magnitude smaller than the contribution due 
to pion exchange Pomeron exchange also does not give a significant contribution since 
neutron production due to diffractive dissociation is believed to be « 6% of the pion exchange 
contribution P]. The present data sample has been used to estimate a 2% diffractive dissociation 
contribution to leading neutron production by determining the fraction of events with a large 
rapidity gap extending into the LAr calorimeter. We have neglected additional backgrounds 
such as neutron production due to resonance decays. 

The pion, pomeron and reggeon flux factors have been determined using hadron-hadron 
data. The pion flux factor f n / p which we have used for neutron production is the same as the 
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one used in reference [[5p: 



fir/p{z,t) = C 



167T 2 



l-z 



exp 



(11) 



where gl Np /{4n) = 13.6 ± 0.3 [0, a' n = 1 GeV~ 2 , i? 2 = 0.3 GeV 2 and the square of 
the Clebsch-Gordon coefficient is C = 2/3. For proton production via n° exchange we use 
the same flux factor with C = 1/3. The pomeron and reggeon flux factors are parameterized 
astHEl]: 



fv/ P (z,t) 



54.4 GeV~ 



'1 _ z y-2a F (t) 



exp 



\2R%t 



(12) 



/«/,(*, t) = (1 - z) 1 -**® exp \2R^\ , (13) 

where a P (i) = (1.08 + 0.25GeV" 2 t) and a n {t) = (0.5 + 0.9GeV" 2 t). The slopes are 
i?p = 1.9 GeV -2 and R^ = 2 GeV -2 respectively. The modulus squared of the reggeon 
signature factor PL which is approximately equal to two, has been absorbed into the reggeon 
coupling and we have not included reggeon-pomeron interference terms in the model. 

The evaluation of the pion flux factor is not without some theoretical uncertainty. It has been 
pointed out that absorptive corrections, generated by double reggeon pion-pomeron exchanges, 
might play an important role in hadronic reactions in contrast to DIS. Since the pion flux factor 
which we have used was determined using pp — > nX data, it might underestimate the flux of 
pions in the proton for DIS reactions by up to w 30% [|46|,|4"7|]. 

The structure functions for the exchanged particles are basically unknown in the low (5 
region and one has to rely on theoretical models. For the pion structure function F£ we took the 
leading order parametrization by Gliick, Reya and Vogt For the reggeon and pomeron 

structure functions we assume Ff 1 = F£ and Ff = (0.026/0. 12)Ff l following the arguments 
given in reference JTOj] . Measurements of the diffractive structure function F^ 3 - 1 only 
probe the pomeron at high (3 {(3 > 0.04) and it is not possible to use these data to fix the 
pomeron structure function Ff at the low (3 values of the semi-inclusive data (/3 < 3 • 10 -3 ). In 
the small region of overlap however, the QCD fits to [@] are consistent with the pomeron 
model used in this paper as will be discussed below. 

The model gives an acceptable description of the neutron and proton data with 0.7 < z < 
0.9, in view of the fact that all particle fluxes and structure functions were taken from the 
literature and that no adjustment was made. The rate of leading neutron production can be 
described entirely by n + exchange. However, proton production requires contributions from 
both f 2 and n° exchange which are roughly in the ratio 2 : 1 from the model. 

'in reference jio[], equation 7 is missing the reggeon signature factor which is given in equation 5 of refer- 
ence . The two publications also use different values of R^. In reference — 1.2 GeV -2 which leads 
to a 12% difference in the values of the p^-integrated reggeon flux factor at z = 0.8. 
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The shaded-band in Figure || shows the prediction for F 2 P( - 3 ^ in which we have replaced the 
pomeron component in our Regge model with the pomeron component determined using the 
QCD fit to F^ 3 - 1 [Q]. In the QCD fit, the pomeron structure function is parameterized at a low 



scale and evolved to larger Q 2 using the leading order DGLAP [HSp equations. The hard-gluon 



leading order result which we have used for F^ (fit 3 in reference [Q]) is only shown in the region 
in which it is valid (3 < Q 2 < 75 GeV 2 and 0.04 < (3 < 1.0) and it has been interpolated from 
\t\ < 1 GeV 2 to pt < 200 MeV in order to allow comparison with the leading proton data. The 
width of the band reflects the uncertainty in the interpolation to the different kinematic region. 
The pomeron flux factor used in reference [|3|] has been evaluated using «p(t) = «p(0) + a' F t, 
where a P (0) = 1.203 ± 0.020 (stat.) ± 0.013 (syst.) and a P = (0.26 ± 0.26) GeV 2 @]. 
The ZEUS measurement of the slope parameter b = (7.2 ±1.1 (stat.)^'g (syst.)) GeV 2 P^], 
where b = 2R$> — 2«p In (1 — z), has also been used. This comparison demonstrates that the 
HI measurements of F 2 ' and F 2 v ' can both be described by Regge phenomenology. 

We use the measurement of F 2 N ^ at z = 0.7 and the integral of the pion flux factor to 
estimate the pion structure function at low Bjorken-x. Assuming that our Regge model of 
leading neutron production is valid, the quantity F 2 ^/T^ can be interpreted as being equal to 
the structure function of the pion where: 

/•imin 

r 7r (z = 0.7)= / U lp {z = 0.7, t)dt = 0.131. (14) 

Figure^ shows F^^/T^ as a function of (3 for fixed values of Q 2 . The data are compared to 
predictions of several parametrizations of the pion structure function [ 09| , |4O| ,pO|-p2|]. The latter 



are only shown in the Q 2 regions in which they are valid. The data are in good agreement with 
the expectations of the GRV leading order parametrization of the pion structure function. 

The quark and gluon distributions of the pion have previously been constrained in the x ^ 0. 1 
region using Drell-Yan data and direct photon production data obtained by up scattering exper- 



iments (see for example Q53|-p7]]). Our determination using F 2 N(3 ^ is the first result which 
constrains the pion structure function at values of x which are more than an order of magnitude 
smaller. Background contributions and possible absorptive corrections [^,^,^7|], which have 



not been taken into account, are expected to only affect the absolute normalization of our result 
since all of the data are at z — 0.7. 



7 Summary and Conclusions 

The semi-inclusive cross sections ep — > epX and ep — > enX have been measured in the 
kinematic region 2 < Q 2 < 50 GeV 2 , 6 • 10~ 5 < x < 6 • 10~ 3 and p T < 200 MeV. Comparison 
of the proton and neutron data in the same kinematic domain shows that the production cross 
section for leading protons is larger than it is for leading neutrons. This result demonstrates that 
leading proton production cannot be entirely due to pion exchange. 

The LEPTO Monte Carlo program, which is based upon soft colour interactions and a string 
fragmentation model, describes the magnitude and the general shape of the neutron data with 
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z > 0.3. It fails however to describe the rate of leading proton production and the rise in 
the semi-inclusive structure function F^ 3 - 1 as a function of Q 2 . The RAPGAP Monte Carlo 
program, which simulates pion exchange, gives an acceptable description of the neutron data 
with z > 0.7 but does not explain the absolute rate of leading proton production. 

The proton and neutron data are equally well described by fits assuming factorization in 
x, Q 2 and z or /3, Q 2 and z. The scaling violations observed in the measured semi-inclusive 
structure functions F^ 3 ^ and F^^ are similar in size and shape and are close to those seen in 
the GRV parametrizations of the inclusive structure functions of the pion and the proton. 

The neutron and proton data are reasonably well described by a Regge model of leading 
baryon production which considers the colour neutral exchanges of pions, pomerons and sec- 
ondary reggeons. The semi-inclusive cross sections for leading neutrons with 0.7 < z < 0.9 
can be described entirely by n + exchange whereas the semi-inclusive cross sections for protons 
with 0.73 < z < 0.88 require n° and f 2 exchange contributions. In our model, the contribution 
due to f 2 exchange is approximately a factor of two greater than the contribution due to it 
exchange. The f3 and Q 2 dependence of the leading neutron data at z = 0.7 are consistent with 
the GRV leading order parton distributions for the pion. 
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Table 1 : The measured values of F 2 (x, Q 2 , z) for protons with p T < 200 MeV. There is an 
additional normalization uncertainty of 5.6% not included in the systematic error. 
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Table 2: The measured values of F 2 '(x, Q 2 , z) for neutrons with p T < 200 MeV. There is 
an additional normalization uncertainty of 5.7% not included in the systematic error. 
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Figure 1: Top: perspective view of the scintillating fibre hodoscopes and the trigger scintillators 
inside a Roman Pot. Bottom: horizontal cross section through one FPS station. 
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Figure 2: Configuration of the HI forward neutron calorimeter. The calorimeter consists of 
interleaved layers of lead and scintillating fibres. A hexagonal module, see inset, is defined 
by coupling 1141 scintillating fibres to a common photomultiplier located at the rear of the 
detector. 
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Figure 3: (a) The observed neutron energy spectrum in proton beam-gas interactions compared 
to the results of a pp — > nX Monte Carlo simulation based upon pion exchange. The Monte 
Carlo simulates the acceptance and response of the FNC. (b) The same proton beam-gas energy 
spectrum compared to the neutron energy distribution observed in DIS interactions. The proton 
beam-gas energy spectrum has not been corrected for the trigger efficiency which is less than 
100% below 300 GeV. All distributions are normalized to the number of events with E' > 
500 GeV. 
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Figure 4: (a) The observed proton energy spectrum compared to a simulation based upon the 
RAPGAP Monte Carlo generator. In (b) and (c) the observed neutron z and p T spectra, inte- 
grated over the entire kinematic range in x and Q 2 , are compared to reweighted Monte Carlo 
data which result from the unfolding procedure used to correct the data for acceptance and mi- 
gration effects. The Monte Carlo distributions are normalized to the total number of events in 
the data. 
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Figure 5: Measurement of F 2 1 ', for protons with p T < 200 MeV, compared to the predictions 
of the LEPTO and RAPGAP Monte Carlo models calculated using GRV leading order parton 
distributions for the proton and the pion respectively. The inner error bars show the statistical 
errors and the full error bars show the statistical and systematic errors added in quadrature. 
There is an additional 5.6% overall normalization uncertainty for the data points which has not 
been included in the full error bars. 
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Figure 6: The semi-inclusive structure function F 2 ( ', for neutrons with p T < 200 MeV, 
compared to the predictions of the LEPTO and RAPGAP Monte Carlo models calculated using 
GRV leading order parton distributions for the proton and the pion respectively. The inner error 
bars show the statistical errors and the full error bars show the statistical and systematic errors 
added in quadrature. There is an additional 5.7% overall normalization uncertainty for the data 
points which has not been included in the full error bars. 
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Figure 7: Normalized scaling violations observed in the proton and neutron data, computed 
at Q 2 = 10 GeV 2 , compared to the expectations derived from GRV parametrizations of the 
inclusive structure functions for the pion and the proton. The scaling violations of the pion 
structure function have been evaluated as a function of (3 (lower scale), whereas for the proton 
structure function they have been evaluated as a function of x (upper scale). 
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Figure 8: The measured values of F 2 { ' and F 2 1 ; with z > 0.7 compared to a Regge model 
of baryon production. The different contributions are labeled for the figure in the inset. The 
neutron data are described by 7r + exchange whereas the proton data are compared to the sum of 
7T°, pomeron and secondary reggeon (/ 2 ) exchanges. The n° contribution, which is not shown, 
is exactly half the 7r + contribution. The shaded-band is explained in the text. 
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Figure 9: F2 W /T n at z = 0.7 plotted as a function of /3 for fixed values of Q 2 . The quantity 
rV is the pt integrated pion flux factor. Within the framework of the Regge model, /T w 
can be interpreted as being equal to the pion structure function FJ. The data are compared to 
different parametrizations of F£ which are only shown in the Q 2 regions in which they are valid. 
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